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Abstract. We explore the feasibility of applying the Christensen-Dalsgaard diagram to real asteroseismic data and provide 
quantitative measures of the uncertainty associated with the results. We also propose a new kind of seismic diagram, based on 
the determination of the locations of sharp acoustic features inside a star. We show that by combining the information about the 
position of the base of the convective envelope or the He ii ionisation zone with a measure of the average large separation, it is 
possible to constrain the unknown chemical composition or the various parameters characterising the physical processes in the 
stellar interior. We demonstrate the application of this technique to the analysis of mock data for a CoRoT target star. 
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1. Introduction 

Seismic observations of distant stars, either from the ground 
or space, is likely to permit the detection of only low degree 
modes of oscillations (f = 0, 1, 2, 3). In this context, it is neces- 
sary to develop or adapt the existing seismic tools of analysis to 
extract the maximum amount of information contained in these 
modes. One approach is through the construction of suitable as- 
teroseismic diagrams - the technique of combining useful fea- 
tures of the frequency spectrum of stars to extract information 
about the physical processes in the stellar interior that affect the 
frequencies. The basic idea behind such diagrams is that they 
are to be constructed in terms of asteroseismic observables and 
must reveal some structural or evolutionary feature of a star 
through a calibration of theoretical stellar models. 

The most well-known asteroseismic diagram is the so- 
called Christensen-Dalsgaard diagram, or the C-D diagram, in 
short (Christensen-Dalsgaard 1988). This diagram exploits the 
fact that while the average large separation of radial modes, 
(Avo) = (Avo(n)) - (v(n,0) - v(n - 1,0)) reflects the gross 
properties of a star, like its mass and radius, the average small 
separation, {do2) = {do2(n)) = (v(n, 0) - v(n - 1,2)) is more 
sensitive to the innermost layers of the star, and therefore to 
its evolutionary state. The averages are computed over several 
high-order modes (the so-called asymptotic range). By deter- 
mining the position of a star on this diagram, through mea- 
surement of the large and small separations of its oscillation 
frequencies, it is possible to determine its mass and age. 

In this work, we study the feasibility of using the C-D dia- 
gram for real asteroseismic data containing finite errors in the 
frequencies. We investigate the accuracy to which the mass and 
age of a star with unknown characteristics (such as the chem- 
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ical composition or the extent of convective overshoot) can be 
predicted from the C-D diagram (Section 2). 

The rest of the paper introduces a new kind of asteroseis- 
mic diagram. The acoustic location of sharp features inside a 
star such as the boundaries of convective zones or ionisation 
zones can be determined using the oscillatory signal in the fre- 
quencies that they produce (e.g., Monteiro et al. 1994; Basu et 
al. 1994; Roxburgh & Vorontsov 1994). Considering that the 
location of these regions cannot be independent of the general 
stratification inside a star, in Section 3 we propose a new di- 
agram, connecting these independent seismic observables and 
the mean large separation. An application of this seismic di- 
agram to the analysis of synthetic data for the star HD49933 
illustrates how it could help substantially in constructing a re- 
Uable model of the star (Section 4). 

This work utilises several grids of stellar models with 
different stellar parameters which were constructed using 
the CESAM evolutionary code (Morel 1997). These models 
used the OPAL equation of state (Rogers & Nayfonov 2002) 
and OR4L opacity tables (Iglesias & Rogers 1996), comple- 
mented by the low-temperature opacity tables of Alexander & 
Ferguson (1994). Convection was described by the standard 
mixing length theory (Henyey et al. 1965) and nuclear reac- 
tion rates were obtained from the NACRE compilation (Angulo 
et al. 1999). Diffusion of helium and heavy elements was not 
considered for the present study. While it is generally accepted 
that some form of diffusive phenomenon is very likely to ex- 
ist in stellar interiors, we believe that its inclusion would not 
alter the results of this work in a qualitative way. The effect 
of incorporating diffusion would be akin to exploring another 
dimension in stellar parameter space - the various calibration 
curves might change in absolute terms, but that would not affect 
the general conclusions drawn from them. 
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The p-mode oscillation frequencies of low degree modes 
{i - 0,1,2,3) were computed for each model under the 
adiabatic approximation, using the Aarhus pulsation pack- 
age, ADIPLS (Christensen-Dalsgaard & Berthomieu 1991). 
Throughout the work, we have used these low degree modes 
which are the only ones relevant for real asteroseismic data. 

2. Sensitivity of the C-D diagram to stellar 
parameters 

The general idea of the C-D diagram has already been pre- 
sented in Section 1. Recently, Roxburgh & Vorontsov (2003) 
have demonstrated that while the frequency separations them- 
selves are affected by the effects of the surface layers, the ra- 
tio of the small separation to the large separation, defined as 
(ro2) = {roiiri)) = {^^), can exclude such effects to a great 
extent. Similar C-D diagrams can be constructed using these 
ratios, instead of the small separations themselves (Floranes et 
al. 2005). A comparison of the classical C-D diagram and its 
variation for ^ = 0, 1 modes (Mazumdar & Roxburgh 2003), 
with small separations as well as the ratio of separations, can 
be found in the online appendix. 

To investigate whether the C-D diagram can indeed be used 
for real asteroseismic data, we compare the sensitivity of the 
tracks on the C-D diagram toward different stellar parameters 
to the uncertainty associated with placing an observed star on 
this diagram. We estimate the errors in (Avq) and (^02) through 
a Monte Carlo simulation where random errors with a specified 
standard deviation (o-y = 10""* v) are added to the exact frequen- 
cies before computing the separations. Such an error margin in 
the frequencies is consistent with the expectations from the cur- 
rent and future asteroseismic space missions, MOST (Walker et 
al. 2003) and CoRoT (BagUn 2003), as well as recent ground- 
based seismic observations (e.g.. Bedding & Kjeldsen 2003). 
The variance in the average values of the large and small sep- 
arations are estimated to be the uncertainties in these quanti- 
ties for a real data set with comparable errors in frequencies. 
We find that the error in determining the average large sepa- 
ration is only ~ 0.02 juHz. Thus, although the individual fre- 
quency values might have errors up to ~ 0.3 /iHz, the average 
value, (Avo) over a suitably chosen range of radial orders is 
quite a robust quantity. On the other hand, the uncertainty in 
the average small separation, (<io2) is ~ 0.0015 //Hz. Very sim- 
ilar estimates for these uncertainties were found by Audard et 
al. (1995). However, the errors in determining these quantities 
would be affected by the total number of modes observed as 
well as the mass and age of the star. At late stages of evolu- 
tion, particularly in more massive stars, the presence of mixed 
non-radial modes reduces the number of modes available for 
such asymptotic analysis. Typically, this causes an increase in 
the error on the average small separations. The error estimates 
quoted here, and shown in Figs. 1 to 3 are, therefore, only rep- 
resentative. 

The theoretical C-D diagram is typically based on stel- 
lar models constructed with a particular choice of parameters, 
principal among which are the mixing length a, the extent of 
convective core overshoot, dov (both a and dov being measured 
in terms of local pressure scale height) and the initial chemical 
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Fig. 1. The variation of the C-D diagram using the average small sep- 
arations, <<io2) (top panel), and the average ratios, (ro2> (bottom panel) 
with the choice of the mixing length parameter, a, is shown. Two sets 
of models have been used - one with or = 1.8 (solid curves) and an- 
other with a = 1.6 (dashed curves) to illustrate the shift in the tracks. 
All other stellar parameters are identical in the two sets of models. 
Representative errorbars corresponding to relative frequency errors of 
10"'* are shown. 

composition (characterised by (Xo,Zo)). In order to investigate 
the effect of these parameters on the C-D diagram, we con- 
structed separate grids of stellar models with different values 
for these parameters. For each of these sets of models, only 
one parameter (e.g., a) is varied at a time, keeping the others 
(viz., doY, Xq, Zq) identical, so that each effect can be separately 
understood. The frequency range for determining the asymp- 
totic average values of the separations is chosen in terms of the 
scaled frequency, v/ ^/(M/MQ)/(R/R^ , where M and R are 
the mass and radius for each model respectively, to be 1.5 mHz 
to 3.5 mHz. This typically corresponds to averaging over 15 
radial orders. 

In Fig. 1, we show the variation in the C-D diagrams of 

{do2) and (ro2) with the mixing length parameter, a. Each curve 
on this diagram (and Figs. 2 and 3) is an evolutionary track 
of a given mass, indicated in solar units at the ZAMS end of 
the track. We consider two sets of models, corresponding to 
a = 1.8 and 1.6. Since the outer convective envelope grows 
thinner for higher masses, naturally, the effect of changing the 
mixing length reduces with mass. For low mass stars, the effect 
is comparable to the errors in {doi) or (roi), but the tracks are 
still close enough to allow an estimate of the stellar mass within 
5%. 

The effect of convective core overshoot on the C-D diagram 
is more dramatic. We show this effect in Fig. 2 by comparing 
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Fig. 2. The variation of tlie C-D diagram using tlie average small sep- 
arations, <do2) {top panel), and the average ratios, (^02) {bottom panel) 
with convective core overshoot is shown. Two sets of models have 
been used - one with <iov = (solid curves) and another with dov = 0.2 
(dashed curves) to illustrate the shift in the tracks. All other stellar 
parameters are identical in the two sets of models. Representative er- 
rorbars corresponding to relative frequency errors of 10 are shown. 




(Ai^o) (mHz) 

Fig. 3. The variation of the C-D diagram using the average small sep- 
arations, (^02) (top panel), and the average ratios, <ro2> (bottom panel) 
with initial chemical composition is shown. Four sets of models with 
different combinations of (Xo,Zo) have been used - (0.70, 0.02) (solid 
curves), (0.73,0.01) (small dashed curves), (0.70,0.03) (long dashed 
curx'es), and (0.73, 0.02) (dotted cwves) to illustrate the shift in the 
tracks. All other stellar parameters are identical in the two sets of 
models. Representative errorbars corresponding to relative frequency 
errors of 10"* are shown. 



models without overshoot (day - 0) and those with Jqv = 0.2. 
The effect of overshoot in increasing the total radius of the star 
at a given evolutionary stage (in terms of the central hydrogen 
abundance, Xa) is clearly seen in the extension of the tracks on 
the C-D diagram to lower values of (Avo). Further, the value of 
(c/oa) also reduces with overshoot, reflecting the higher age of 
the star. In effect, the tracks on the C-D diagram deviate consid- 
erably from the standard ones when overshooting is included 
in the models. It is evident, therefore, that we need careful cal- 
ibration of the effect of overshoot if one is to estimate the mass 
and age of a star (with mass ^ LIM©) directiy from the C-D 
diagram. 

While the metallicity of a target star is often known from 
spectroscopic data, its helium content is typically unknown, 
thereby one is left with the ambiguity of one unknown param- 
eter concerning the chemical composition. We consider three 
different scenarios which approximately reflect the various pos- 
sibilities for the chemical composition of a star with a known 
metallicity (or [Fe/H]). We compare the tracks on the C-D di- 
agram for all these cases with respect to those for the standard 
solar composition in Fig. 3. 

The first case ((Xq.Zq) = (0.73,0.01)) corresponds to a 
typical metal-poor star ([Fe/H] ~ -0.25) with proportionally 
reduced helium content. This would exemplify the situation 
where we would only have an indirect estimate of the helium 



abundance from the metallicity, based on assumptions about 
the chemical evolution of the galaxy. The second combina- 
tion {{Xq, Zo) - (0.70, 0.03)) mimics a star which is metal-rich 
([Fe/H] ~ 0.25), but has near-solar helium abundance. The 
third case ((Xo.Zq) = (0.73,0.02)) reflects the chemical com- 
position of a star which has solar metallicity, but is known to 
be over-abundant in hydrogen. It is evident that in each case, 
the tracks depend substantially on the chemical composition. It 
appears that even a slight change in the balance of helium and 
heavier elements at a constant hydrogen abundance produces a 
large shift in the C-D diagram. In summary, the effectiveness of 
the C-D diagram may be strongly undermined if the chemical 
composition is not fully known. 

The importance of the equation of state (EOS) in the deter- 
mination of internal stellar structure and its seismic manifes- 
tation cannot be ignored. However, by comparing two sets of 
models with OPAL and EFF (Eggleton et al. 1973) EOS we find 
that the shift in the tracks are at a level much lower than that 
induced by the other parameters discussed above, especially at 
higher masses. Therefore, it would be fair to deduce that the 
choice of the EOS would not greatly affect the extraction of the 
mass and age of a star from the C-D diagram. 
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3. A-t diagram 

3. 1. Probing sharp acoustic features 

While the C-D diagram is indeed a powerful tool to estimate 
the mass and age of a star from its oscillation frequencies, the 
process of averaging the frequency separations actually reduces 
the information content of the frequency spectrum. One of the 
other means to investigate the structure of the stellar interior 
from the oscillation frequencies without resorting to explicit 
modelling of the star is to utiUse the oscillatory signal in the 
frequencies to determine the acoustic depth 




of a sharp feature, lying at a radius of r^, c being the sound 
speed and R the total radius of the star (e.g., Gough 1990). 

The acoustic depth of the He n ionisation zone (hereafter 
referred to as HIZ) and the base of the convective envelope 
(hereafter referred to as BCZ) can be easily estimated using 
the second differences of the frequencies which are defined 
as 6^Vn,e = Vn+i/ - 2v„/ + v„_i/ (Mazumdar & Antia 2001). 
Following Basu et al. (2004), we fit a function of the form 

6^v = ^ai +a2V+ 

+ ^^1 + ^ + ^ j sin(4;rvfHiz + (f>mz) 

(C2 C3 \ 

Cl + — + ^ j sin(4;ryfBcz + 0bcz), (2) 

to the second differences of the frequencies. The parameters 
ai,bi,Ci (i = 1,2,3) and thiz, tbcz, <^hiz and ^bcz are ob- 
tained from a least-squares fit. The frequency range chosen for 
performing the fit is 1.5-3.0 mHz, after scaling by the factor 
(Ayo)/(Avo)o, where (Ayo) and (Ayo)© (- 135/iHz) are the av- 
erage large separations for the concerned star and the Sun re- 
spectively. 

However, it is generally found that the acoustic depths of 
the BCZ and HIZ derived from the oscillatory signal in the fre- 
quencies (tbcz and thiz) are systematically higher than their 
true values (tbcz and thiz) calculated directly from the model 
according to the definition, Eq. (1). This is due to the fact that 
the estimation of the acoustic depth of the layer of disconti- 
nuity through the oscillatory signal in frequencies is biased by 
the surface phase shift (Vorontsov et al. 1991). However, as 
pointed out by Christensen-Dalsgaard et al. (1995) and Ballot 
et al. (2004), this bias can be largely eliminated in the estima- 
tion of the acoustic radii, fscz and fniz of the BCZ and HIZ, 
instead of the acoustic depths: 

^Bcz = S) - Tbcz . and fniz = ?o - thiz . (3) 
where fo = 2<av„> is an estimate of the total acoustic radius of 

the star, fo = Jjf drjc. 

By estimating the acoustic radii, Fbcz and fHiz> for several 
models of different mass and age by this method, we find that 
the true acoustic radii and the estimated values differ by ^ 5% 
- a remnant effect of the surface phase shift. The differences 
tend to increase with mass, because for stars more massive than 



~ 1.5Mo the outer convective envelope becomes thin and the 
BCZ and the HIZ lie close to each other. In such a situation, 
it is difficult to fit the two oscillatory components separately. 
Therefore, this method would not work very well for stars more 
massive than I.5M0. 

This difference between the model values of the acoustic 
radii and their estimated values from the oscillatory signal ex- 
ists even for exact frequencies and might be considered as a 
systematic bias in the process. In addition, there would be a 
certain random error in the case of actual frequency data asso- 
ciated with observational uncertainty. We have estimated this 
error through Monte Carlo simulations similar to that described 
in Section 2. We find that the uncertainty in the estimates of 
the mean density-scaled acoustic radii of the BCZ and the HIZ 
is typically around 40 s for relative frequency errors of 10 
Since the uncertainty in the average large separation turns out 
to be tiny (Section 2), practically the entire contribution to- 
wards this error is from the fitting of Tbcz and thiz- This is 
similar to the error margins provided by Basu et al. (2004), 
and comparable to that for the BCZ estimated by Ballot et al. 
(2004). 

3.2. A new seismic diagram 

We have seen in Section 2 that the average large separation, 
which is an index of the gross properties such as the mean den- 
sity or the total acoustic size, can be determined to a fairly good 
accuracy from an observed frequency spectrum of a solar-type 
star. Further, in Section 3 . 1 , we established that the acoustic po- 
sition of the BCZ and the HIZ can be determined to an accuracy 
of at least 95 % from the oscillatory signal in the frequencies de- 
spite the presence of some systematic bias due to unknown sur- 
face effects. Importantly, these estimates are obtained without 
explicit modelling of the stellar structure. However, the acous- 
tic location of these layers inside the star cannot be indepen- 
dent of the general stratification or the chemical composition 
of the star. Therefore, the model-independent estimates of fecz 
and tmz indeed provide clues about the internal structure of the 
star. 

In order to combine these independent pieces of informa- 
tion about the gross properties and the detailed structure of a 
star, we construct a diagram connecting the average large sep- 
aration and the acoustic radius of the BCZ or the HIZ. The ex- 
ample of such a diagram, which we shall call the A-t diagram, 
is given in Fig. 4. Each Une on this diagram is actually a locus 
of the acoustic radius of the BCZ (or the HIZ) as a function of 
the average large separation as a star of a given mass evolves on 
the main sequence. Thus it may be considered as a kind of seis- 
mic diagram, analogous to the C-D diagram. In this diagram, 
the ordinates are actually not the true acoustic radii Ibcz and 
tmz, but rather are the values derived from the acoustic depths 
Tbcz and thiz (model values) and the mean large separation 
(Ayo) (derived from exact frequencies), through Eq. (3). This 
is done deliberately to minimise the systematic error described 
in Section 3.1 in the calibration process. A diagram similar to 
that of Fig. 4 may be constructed using the acoustic radii turn, 
instead of fncz- Since the properties of these diagrams are, in 
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Fig. 4. An example of the A-t diagram, wiiere tlie acoustic radius of 
the BCZ is plotted as a function of the average large separation. The 
ordinate is plotted in logarithmic values for the sake of clarity. Each 
curve on this diagram is an evolutionary track of the indicated mass in 
solar units. 



general, quite similar, it is enough to focus on the A-Fbcz in the 
following discussion. 

3.3. Sensitivity of A-t diagram to stellar parameters 

It is evident from Fig. 4 that the A-f diagram has the diagnos- 
tic capabihty to predict the mass and evolutionary status of a 
star, in a manner similar to that of the C-D diagram. However, 
the tracks on this theoretical diagram would naturally depend 
on all the parameters involved in the stellar modelling. Direct 
determination of the mass and age would only be possible if 
the other stellar parameters are known a priori. Since typically 
this is not the case for real stars, we have investigated the sen- 
sitivity of the A-t diagram to stellar parameters by using sev- 
eral grids of models, for four different masses, 0.8, 1.0, 1.2 and 
1.5Mq. Since the acoustic radii Ibcz and ?hiz cannot be reliably 
obtained for more massive stars, they are not relevant for this 
exercise. 

In Fig. 5, we show the A-t tracks for models with two dif- 
ferent values of the mixing length parameter, a; = 1.8 and 1.6. 
The errorbars for t^cz for relative frequency errors of 10""* are 
magnified by a factor of 3 in these diagrams to increase visi- 
biUty. We find that the tracks for tscz as a function of (Avq) 
are shifted almost parallelly as a is changed. The amount of 
this shift depends on the mass - while at lower masses it is too 
small to be detected given the uncertainty in the seismic deter- 
mination of fficz, at masses ^ \.2Mq, the shift is quite large. 
This would imply that this diagram is sensitive to a at higher 
masses, and might provide at least an indicative value of the 
same, provided other parameters are known. 

The effect of convective core overshoot is expected to be 
small on these A-t diagrams, since the BCZ Ues far apart from 
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Fig. 5. The sensitivity of the A-f diagram with respect to the mixing 
length parameter a is illustrated. For each mass, indicated in the dif- 
ferent panels, two sets of models have been used - one with or = 1.8 
{solid curves) and another with a = 1.6 {dashed curves). AU other 
stellar parameters remain identical. Representative 3cr errorbars for 
the values of the acoustic radii for a relative error of lO""* in frequen- 
cies are shown. 
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All other stellar parameters remain identical. Representative 3cr error- 
bars for the values of the acoustic radii for a relative error of 10"'' in 
frequencies are shown. 



the convective core. Figure 6, however, confirms this only par- 
tially. We find that fecz is indeed sensitive to c/ov at later stages 
of evolution for masses between 1.2 and \ .5Mq. 

Figure 7 illustrates the effects of chemical composition on 
the A-t diagram. The compositions chosen for comparison to 
the standard one are identical to the ones described in Section 2. 
We find that the BCZ location is sensitive to the chemical con- 
tent in all three scenarios. The effect is very pronounced at 
higher masses, which would mean that the A-t diagram may 
be used to constrain the chemical abundances, at least for 1.2- 
I.SMq stars. 
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Fig. 7. The sensitivity of tlie A-t diagram to initial chemical com- 
position is illustrated. Four sets of models with different combi- 
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identical. Representative 3cr errorbars for the values of the acoustic 
radii for a relative error of 10 in frequencies are shown. 



4. Application of the A-t diagram: The hare & 
hound exercise for HD49933 

There might be two approaches to using the A-t diagram. The 
first is to determine any one parameter among M, (Xq, Zq), a 
or doy directly. But this would require independent knowledge 
about the remaining parameters, since each of these parameters 
affect the A-t diagram considerably. Some constraints about 
the mass, age and the chemical composition might be indepen- 
dently gathered in case the target star is a member of a binary 
system or a cluster. In such cases, it might be possible to use 
the A-t diagram to estimate the mixing length or the overshoot 
parameter. 

However, we can also envisage the use of the A-t diagram 
in another way, when none of the stellar parameters are known 
a priori. In this case, given an initial trial model, one can test 
its relative position on the A-t diagram with respect to the 
observed data to gain insight into how a particular parameter 
needs to be changed from the initial guess in order to find a 
better match with the data. This process has to be iterative, 
but would considerably reduce the task of searching in a multi- 
dimensional parameter space for the best model. 

In this section, we illustrate the use of the A-t diagram 
as applied successfully in the seismic modelling of a CoRoT 
primary target star, HD49933, based on synthetic data in the 
course of a hare and hound exercise. This exercise, carried 
out within the CoRoT Seismology Working Group, was con- 
ducted in four principal steps - (i) construction of a the- 
oretical model for a star, based on its position in the HR 



diagram and available information about its chemical com- 
position, and calculation of its low-degree eigenfrequencies, 
(ii) conversion of the frequencies into a time series, taking into 
account the duration of CoRoT observation (150 days) and the 
expected instrumental and background stellar noise, (iii) ex- 
traction and mode identification of frequencies from this simu- 
lated time series, without knowledge of the original exact fre- 
quencies, and (iv) seismic interpretation of the extracted fre- 
quencies without knowledge of the original input model. Each 
of these steps were carried out by independent groups, main- 
taining confidentiality of the original input data, thereby sim- 
ulating a situation which would be encountered for real data 
gathered from the CoRoT mission. The general description 
of this exercise can be found on the web (http://www.ias.u- 
psud.fr/virgo/html/corot/datagroup/hh.html) and here we shall 
restrict ourselves only to the particular case of the seismic mod- 
elhng of the star HD49933i. 

The basic parameters available for HD49933 are sum- 
marised in Table 1 (stage 0). The data on which the seismic 
interpretation was based consisted of 67 individual frequencies 
of degree { = and 1, lying between 100 and 3100 juHz, ex- 
tracted in step (iii) described above. The errors in the frequen- 
cies ranged between 0.02 and 0.46 juHz, with the highest errors 
reported for intermediate frequencies (~ 2000 /iHz). The aver- 
age large separation, (Avq), was found to be 90.4 ± 0.2 //Hz. 
Since the mean density of the star is not known at this stage, 
the average large separation was calculated for several differ- 
ent limits in the entire asymptotic range, and the variation in 
its value was included as an additional error. The mean den- 
sity was estimated to be 0.465 ± 0.015 pg through a calibration 
of the average large separation. Unfortunately, only a few fre- 
quencies for { - 2 and 3 could be extracted from the simulated 
time series, and therefore, we had no small separations doi to 
utilise in our seismic analysis. This automatically ruled out the 
use of the C-D diagram. 

We derived the acoustic depths tbcz and thiz by fitting an 
expression, similar to that of Eq. (2) to the second differences 
of the frequencies, as shown in Fig. 8. For this exercise, we 
used all modes between 700 and 2700 juHz, the range being 
chosen to use as many modes as possible, while not deviat- 
ing too much from the standard choice as explained by Basu 
et al. (2004). This yielded values of t^cz = 4085 ± 68 s and 
tmz = 4866 ± 47 s. The comparison of these results to the orig- 
inal model values are shown in Table 1 (stage 1). The excellent 
agreement of these values, obtained prior to any modelling of 
the star, indeed confirms the diagnostic power of the oscillatory 
signal in the frequencies. In practice, however, the original val- 
ues were not revealed at this stage, but only at the end of the 
entire exercise. 

In the next step, we proceeded to use the A-f diagrams to 
arrive at the closest model for the star, already taking into ac- 
count the constraints imposed on the stellar parameters by the 
position of HD49933 on the HR diagram. In the top panel of 



' The relevant work in steps (i), (ii) and (iii) for HD49933 was car- 
ried out by 1. Roxburgh, C. Barban and T. Toutain respectively. The 
seismic interpretation in step (iv) described here was carried out by 
the author in collaboration with E. Michel. 
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Fig. 8. The oscillatory signal in the second differences of the frequen- 
cies of HD49933 (shown as data points with respective errorbars) are 
fitted by the expression (2) (solid curve) to extract the values of tbcz 
and thiz. 



Fig. 9, we show the evolutionary tracks of three sets of models 
with the same mass, overshoot and the mixing length parame- 
ter, but differing in chemical composition, and passing through 
the errorbox of HD49933. The A-t diagram is illustrated in the 
bottom panel, where the range in the values o/(Avo) and Fbcz 
for each track are restricted to correspond to their overlap with 
the box on the HR diagram. The box on the A-t diagram repre- 
sents the values of (Avq) and fecz that we have already derived 
from the frequencies. This is an example of one of the A-f dia- 
grams that we utilised in this analysis, in each case varying one 
stellar parameter. 

The significant result to be noted here is that a much 
stronger constraint on the varying parameter (in this case the 
chemical composition) and the age can be obtained by invok- 
ing the A-/ diagram than is possible through the errorbox on 
the HR diagram alone. Of course, this is not the definitive set of 
parameters for the best possible model, since we have allowed 
only the chemical composition to vary, while fixing the other 
parameters. We need to change all the parameters one by one 
within the range allowed by the HR diagram, and test for the 
remaining parameter through the A-t diagram. But the strength 
of the A-; diagram lies in the fact that starting from an initial 
guess model, we can get a very good indication as to in which 
direction we should move in the multi-dimensional parameter 
space to arrive at the best model. The task of searching for the 
best possible combination of parameters is greatly reduced by 
the A-f diagram through the additional information about the 
location of the BCZ and HIZ. 

After several iterations, we were able to converge to a set 
of models that satisfied the constraints on the A-t diagram for 
both fscz and fniz. The final best models were selected on the 
basis of comparisons of large separation values between the 
data and the models. The absence of small separations in the 
data prevented us from constraining the parameters further. We 
estimated errors on each parameter to represent the range of 
values that would produce models which are all consistent with 
the errorbox on the A-t diagram, as well as the individual large 




3.86 3.84 3.82 

log T.,, (in K) 



3.80 




4400 - 



4200 



4000 



3800 



85 90 
(Ai/o) (mHz) 

Fig. 9. Top panel: The position of HD49933 is shown on the HR dia- 
gram. Three evolutionary tracks for models with a given mass, over- 
shoot and mixing length parameter, but different chemical composi- 
tions are illustrated. Bottom panel: The A-f diagram for Fbcz is shown 
with the three tracks only for the portion where they lie inside the box 
in the top panel. The derived values of (Avq) and fscz for HD49933 
are represented by a box. 



separation values. Needless to add, these models are already 
consistent with the hmits imposed on the HR diagram. 

The parameters for our best models are shown in Table 1 
(stage 2), alongwith those used for the original input model. 
Clearly, a remarkably good match with the input model was 
achieved. We should add a caveat though - as it turned out, both 
the input modelling and the seismic analysis had utihsed the 
same evolutionary code, CESAM; this might have contributed 
partially to the exactness of the result. However, the important 
point is not the goodness of the fit, but rather the illustration of 
the use of a new seismic technique, the A-f diagram, which led 
us to the final result. 

5. Conclusions 

We have investigated the possibility of using the C-D diagram 
for real asteroseismic data by comparing the error in the sep- 
arations propagated from the observed frequencies against the 
sensitivity of the tracks toward different stellar parameters. The 
largest uncertainty in caUbrating the C-D diagram arises from 
the convective core overshoot parameter and the chemical com- 
position of the star. However, the effect of the mixing length 
parameter and the equation of state are found to be relatively 



8 



Anwesh Mazumdar: Asteroseismic diagrams for solar-type stars 



Table 1. The input parameters and results from a seismic analysis of 
HD49933 in a mock exercise are shown. The table is divided in three 
parts, classified by stages of the seismic analysis: (0) initial constraints 
provided to both modelling team and seismic interpretation team; (1) 
results obtained from seismic analysis, prior to modelling; (2) final 
results obtained from seismic modelling. 



Stage 


Quantity 


Non-Seismic 


Seismic 


Original 






conslrainl 


value 


input 





Mm 
[Fe/H] 


6700 + 100 K 
3.35 ±0.10 mag 
-0.32 ±0.10 






1 


P/Po 




0.465 ± 0.015 


0.467 




fscz 




4085 ± 68 s 


4031 s 




'hiz 




4866 ± 49 s 


4805 s 


2 


M/Me 




1.24 ±0.01 


1.25 




Xc 




0.495 + 0.010 


0.50 




Xo 




0.72 + 0.01 


0.73 




Zo 




0.01 ± 0.005 


0.01 








0.20 ±0.01 


0.20 




a 




1.60 ±0.05 


1.50 



small. We estimate that if the chemical composition is known 
independently, it would be possible to determine the mass and 
age of a low mass star within 5% using the average frequency 
separations. For a star with a convective core, the uncertainty 
would increase due to the unknown extent of convective core 
overshoot. 

We propose a new kind of seismic diagram (which we call 
the A-f diagram) for low mass main sequence stars, connect- 
ing the acoustic location of sharp features inside a star and its 
mean large separation. Since the acoustic location of the base 
of the convective envelope or the He ii ionisation zone can be 
estimated to a precision of ~ 95%, the position of a star on 
this diagram can be determined fairly accurately. The relation 
of such sharp features with the internal stratification of the star 
and the dependence of the mean large separation on its gross 
properties constitutes the physical basis behind this diagram. 
This diagram is designed to help in modelling the star by ex- 
ploiting the sensitivity of the acoustic location of the sharp fea- 
tures toward different stellar parameters. 

The use of the A-f diagram is illustrated through the seis- 
mic interpretation of simulated data for a CoRoT primary tar- 
get star, HD49933, carried out in the form of a hare and hound 
exercise. Firstly, independent of comparison of actual models, 
the acoustic radii of the base of the convective envelope and 
the He ii ionisation zone are recovered to a good accuracy by 
analysing the oscillatory signal in the second differences of the 
frequencies. Further, all the stellar parameters used in the orig- 
inal input model could be correctly estimated through the seis- 
mic modelling based on the A-f diagram. 

It is found that the A-? diagram can constrain the multi- 
dimensional space of stellar parameters substantially and pro- 
vide an iterative method of quickly converging to the best 
model, starting from an initial guess. The success of the method 
could be somewhat sensitive to the position of the target star on 
the HR diagram, since the effect of the various stellar param- 
eters on the A-t diagram is not uniform with respect to mass 



or age of the star. We find that at least in the mass range of 
~ 1-1 .5Mq, this method is likely to be powerful. The efficiency 
can be further enhanced if the A-t diagram can be used in con- 
junction with the C-D diagram, which was unfortunately not 
possible in this particular hare and hound exercise. 

The application of asteroseismic diagrams to oscillation 
frequencies of solar-type stars looks to be quite plausible, pro- 
vided the dependence of these diagrams on different unknown 
stellar parameters is taken into account in the calibration pro- 
cess. On the other hand, this sensitivity of seismic diagrams Uke 
the A-f diagram can actually be exploited to construct reliable 
models for target stars based on the model-independent results 
obtained from their oscillation frequencies. Clearly, more tests 
are required with real data to refine these methods. 
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Fig. 10. The classical C-D diagram, plotting the average small sep- 
aration, ((io2>> against the average large separation, (Avq), for a grid 
of main sequence models between mass O.SMg and IOMq. The bold 
lines follow the evolution of each stellar model with a given mass, 
which is indicated at the top of the track. The dotted lines connect 
the models at same stages of evolution, indicated by the central hy- 
drogen abundance, Xc. Representative errorbar for (doi) of a IMq star, 
corresponding to relative frequency errors of 10"'* is shown. 



6. The classical C-D Diagram and its variations 

We present a "classical" C-D diagram in Fig. 10, involving the 
large separations (Avq), and the small separations {^02). The 
stellar models used for this diagram were constructed using the 
OPAL equation of state and had an initial chemical composi- 
tion of (Xo,Zo) = (0.70,0.02). The mixing length parameter 
(in terms of the local pressure scale height), a, is 1.8 and no 
convective core overshoot was considered. 

In Fig. 11, the corresponding C-D diagram using the ra- 
tios of separations is shown. As pointed out by Christensen- 
Dalsgaard (2004), the filtering out of the surface effects indeed 
increases the orthogonality of the two effects regarding mass 
and age. Further, we find that this C-D diagram has improved 
"resolution" at higher masses, i.e., the lines corresponding to 
masses greater than 2Mq do not collapse as much as in the 
classical C-D diagram (cf. Fig. 10). This implies that the diag- 
nostic capability of the C-D diagram can both be improved and 
extended to higher masses if we use the ratios of separations 
instead of the separations themselves. 

The kinks seen in the tracks for models with mass 1.2- 
I.SMq are not results of inconsistencies in the models, but are 
in fact related to the evolution of the convective core. In this 
mass range, the convective core grows in mass (or radius) for a 
significant part of the main sequence evolution, before begin- 
ning to shrink later. The kink corresponds to the evolutionary 
phase where this phase of growth of the convective core ends 
and the shrinking phase begins. This effect is even more pro- 
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Fig. 11. A C-D diagram using the average ratio, (roa), of the small 
separation to the large separation. The stellar models used are identical 
to those in Fig. 10. 

nounced in models including convective core overshoot (see 
Fig. 2). 

The construction of the C-D diagram illustrated above re- 
quires the determination of a set of /' = mode frequencies as 
well as ^ = 2 frequencies of similar radial orders. In practice, 
however, we might encounter the situation where only ^ = 
and I - \ modes are identified for a target star, but not the 1-2 
modes. Following Eqs. (4) and (5) of Roxburgh & Vorontsov 
(2003), we can define the 5-point separations between these 
two sets of modes. The C-D diagrams with these separations 
and the corresponding ratios (Eq. (8) of the same paper) are 
shown in Fig. 12 and Fig 13. 

It is interesting to note the difference in the structure of 
the two diagrams in Figs. 10 and 12. The C-D diagram with 
the {ddo\) separations undergoes a "twist" at higher masses. 
This is, in fact, to be expected from the behaviour of these 
separations as functions of frequency. In Fig. 14 we plot both 
do2 and ddoi as functions of the scaled frequency, v/ sjp/po, 
where p denotes the mean density, for a sequence of models 
of mass 1.5 Mo at different evolutionary stages on the main se- 
quence. The range of scaled frequencies used to compute the 
average values in Figs. 10 and 12 are shown as well. We find 
that while the value of do2 decreases monotonically with age 
at all frequencies, at low frequencies the value of ddoi actually 
increases with age, flattening out to an almost constant value at 
intermediate frequencies. At frequencies close to the acoustic 
limit ddo\ decreases slightly with age. This explains the "twist" 
in Fig. 12. Evidently, the range of frequencies over which the 
average values are computed plays a crucial role in determin- 
ing the structure of the C-D diagram. The physical reason for 
this difference of behaviour between do2 and ddoi, as hinted 
by Audard et al. (1995), possibly lies in the difference in the 
position of the iimer turning points of the ^ = 1 and { = 2 
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Fig. 12. An alternative C-D diagram using the average small separa- 
tions {ddoi)- The stellar models used are identical to those in Fig. 10. 
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Fig. 13. A C-D diagram using the average ratio, (roi), of the small 
separation to the large separation. The stellar models used are identical 
to those in Fig. 10. 
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modes close to the stellar core. While the lowest frequency 
( = I modes would be reflecting off a layer close to the peak 
in the Brunt-Vaisala frequency at the boundary of the convec- 
tive core, the { = 2 modes may not be penetrating that far. The 
difference in behaviour of these two smaU separations bears a 
clue to the evolution of the Brunt-Vaisala frequency with age. 



